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By solving self-consistently an effective Hamiltonian including interactions for both antiferromag- 
netic spin-density wave (SDW) and d-wave superconducting (DSC) orderings, a comparison study 
is made for the local magnetic structure around superconducting vortices and unitary impurities. 
To represent the optimally doped regime of cuprates, the parameter values are chosen such that the 
DSC is dominant while the SDW is vanishingly small. We show that when vortices are introduced 
into the superconductor, an oscillating SDW is induced around them. The oscillation period of the 
SDW is microscopically found, consistent with experiments, to be eight lattice constants (8ao). The 
associated charge-density wave (CDW) oscillates with a period of one half (4ao) of the SDW. In the 
case of unitary impurities, we find a SDW modulation with identical periodicity, however without 
an associated CDW. We propose neutron scattering experiments to test this prediction. 
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The common feature of all high-transition-temperature 
(high-T c ) cuprates is the proximity between antifer- 
romagnctic (AF) and d-wave superconducting (DSC) 
phases controlled by the doping — At low temperatures, 
the AF order of the parent compounds is replaced by 
a DSC order upon doping. The main theme of high- 
Tc superconducivity research centers on how to estab- 
lish the connection between these two kinds of order- 
ings, i.e., whether they exclude each other, or coexist 
microscopically. In one of the strategies, the interplay 
between the AF and DSC ordering is explored by weak- 
ening the superconductivity in the optimally or slightly 
ovcrdopcd regime. The inelastic neutron scattering (INS) 
measurements on both optimally doped and underdopcd 
La2_ 2; Sr 2 ;Cu04 (LSCO) samples by Lake and co-workers 
provide the first evidence of the field-induced magnetic 
order |[| ^j. The field- induced spin fluctuations have a 
spatial periodicity of 8ao with the wave vector pointing 
along the Cu-0 bond directions. The principle mag- 
netization oscillations in the vortex state arc coherent 
over a distance Lm > 20ao, substantially longer than 
the superconducting coherence length £o (~ 5ao). The 
field-induced enhancement of the Bragg peak intensity 
was also observed by Khaykovich et al. || in the elas- 
tic neutron scattering (ENS) measurement on a related 
material, La2Cu04+ y (LCO), but with Lm > lOOao, in- 
dicating the existence of the field-induced static AF or- 
der of 8ao periodicity. Strong AF fluctuations have also 
been observed by Mitrovic et al. B in a high-field nuclear 
magnetic resonance (NMR) imaging experiment on near- 
optimally doped Y^^Cn^Oi-x (YBCO). More recently, 
the scanning tunneling microscopy imaging by Hoffman 
et al. has revealed the quasiparticle states around 
the vortex cores in slightly overdoped Bi 2 Sr2CaCu208 +2 , 
(BSCCO) a Cu-0 bond-oriented "checkboard" pattern 
with 4czo periodicity. The periodicity (4ao) of charge 
modulation is one half of that (8ao) of the field- induced 
SDW modulation. Theoretically, Dcmler and co-workers, 



by focusing their attention on the far-field regions outside 
the vortex core, proposed a phenomenological model || 
that when the superconductivity is weakened by the cir- 
culating currents caused by the vortices, a coexisting 
SDW plus DSC phase appears surrounding the core; 
more recently they also argued JjJ that the halved peri- 
odicity of the static CDW modulation is associated with 
the "Friedel oscillation of the spin gap" . As an extension 
of earlier work within the SO(5) theory ||, ^J, originally 
suggesting an AF insulating region inside the core, Hu 
and Zhang showed that the vortex-induced AF re- 
gion can be greater than due to the light effective 
mass of the dynamic AF fluctuations at optimal dop- 
ing. In these two cases, the 8ao periodicity of the mag- 
netic modulation is not fully understood due to the phe- 
nomenological nature of the models. Distinctly, stripe 
models Q ^ ^ Q predict that the spin modulation 
of wavelength A in cuprates should be associated with the 
charge modulation of wavelength A/2. However, two di- 
mensional modulation is apparent in the STM image [gj. 
One of us and C.S. Ting |1JJ have applied an effective 
microscopic mean-field model accounting for the compe- 
tition between the AF and DSC orderings to success- 
fully explain the peak-split structure around the Fermi 
surface in the local density of states (LDOS) at the vor- 
tex core center Jl7|, [l8| [Hj). Nevertheless, due to the use 
of the simplest band structure parameter values, the spin 
and charge modulation is too weak for the determination 
of the periodicity. In this work, using realistic model 
parameter values, we present a comprehensive study of 
the spin/charge structure around the vortices and non- 
magnetic unitary impurities in optimally doped high-T c 
superconductors. We find, for the first time, microscop- 
ically: (i) For the vortex case, the oscillation periods of 
the SDW and the associated CDW are indeed 8ao and 
4ao, respectively, in good agreement with the INS M and 
STM JEj experiments, (ii) Around the nonmagnetic uni- 
tary impurity, the modulation of the SDW order still cx- 



2 



hibits 8ao periodicity. However, the charge density shows 
only the Friedel oscillation with a period of the Fermi 
wavelength due to the strong potential scattering from 
the impurity. 

Consider a minimal model defined on a two dimen- 
sional (2D) lattice square , in which the on-site repulsion 
is solely responsible for the antiferromagnetism while the 
nearest neighbor attraction causes the <i-wave supercon- 
ductivity |ll| . With the application of an external mag- 
netic field and/or in the presence of nonmagnetic impu- 
rities, the effective mean-field Hamiltonian can be diag- 
onalized by solving sclf-consistently the Bogoliubov-dc 
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Here (m™,,^) is the quasiparticle wavefunction corre- 
sponding to the eigenvalue E n , the single particle Hamil- 
tonian Ttij,o- = —tije lipi3 + (mis + €j — n)6ij, where 
tij = t for the nearest neighbor hopping while Uj = t' 
for the next-nearest neighbor hopping, ti is the single 
site potential describing the scattering from impurities, 
and the Peierls phase factor (py = ^ J r * A(r) • dr with 
$0 = hc/2e the superconducting flux quantum in the 
presence of an externally applied magnetic field. Notice 
that the quasiparticle energy is measured with respect to 
the Fermi energy. The self-consistency conditions read: 
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where U, V are the strength of the on-site repulsion and 
the nearest neighbor attraction, respectively, and the 
Fermi distribution function f(E) = l/[e E l kBT + 1]. Here 
the summation is over the cigenstates with both positive 
and negative eigenvalues jl5|. We report results below 
for two cases at zero temperature. For the Abrikosov 
vortex state, the magnetic field effects enter through the 
Peierls phase factor ipij and no impurities arc introduced 
(e, = 0). For the effects of a single impurity, we set A = 
so that ifiij = 0. Hereafter we measure the length in units 
of the lattice constant ao and the energy in units of the 
hopping integral t. To mimic a hole-like Fermi surface, as 
relevant to the hole-doped cuprates, we take t' = —0.2. 
As a model calculation, the on-site repulsion and pairing 
interactions are taken to be, U = 2.5 and V = 1.0. In 
addition, the filling factor rif = a rii a /N x N y is fixed 
to be 0.84, which corresponds to an optimal hole dop- 
ing rih = 1 — rif = 0.16. Here N x ,N y are the linear 



dimensions of the unit cell under consideration. We use 
an exact diagonalization method to solve the BdG equa- 
tion (Q) sclf-consistently. In the absence of magnetic field 
and impurities, we recover all results reported in |T(j] . In 
the optimal doping regime (n/j = 0.16), the AF SDW or- 
der is absent while the DSC order is homogeneous. When 
the DSC is weakened by the application of an external 
magnetic field or by introducing impurities, the AF SDW 
will be nucleated in the region where the DSC order is 
depressed. 

Field-induced SDW and CDW. When an external mag- 
netic field is applied perpendicular to the 2D Cu-0 plane, 
H = Hi (H c i <C H <C H C 2) : an Abrikosov vortex state 
is formed. As the vortex core is approached, the DSC 
order parameter vanishes topologically due to the circu- 
lating supercurrent surrounding the core. It is assumed 
that the superconductor is in the extreme type-II limit 
where the Ginzburg-Landau parameter k = Ao/£o goes 
to infinity so that the screening effect from the supercur- 
rent is negligible. We choose a Landau gauge to write 
the vector potential as A = (—Hy, 0,0), where y is the 
y-component of the position vector r, so that the Peierls 
phase factor ip^j is uniquely determined. By taking the 
strength of magnetic field, H = 2$o/N x N y , such that 
the flux enclosed by each unit cell is twice <f>o, we solve 
self-consistently the BdG equation (Q) with the aid of the 
magnetic Bloch theorem as given by Eq. (3) in Ref . [jl5| . 
For the calculation, we consider the magnetic unit cell 
of size N x x N y = 48 x 24, which gives a square vor- 
tex lattice. The numerics shows that each unit cell ac- 
commodates two superconducting vortices each carrying 
a flux quantum <I>0! which conforms to the above pre- 
scription for the magnetic field strength. The unit cell 
is equally partitioned between these two vortices, each 
located at the center of area ^ x N y sites. Typical 
results on the structure around one vortex core are dis- 
played in Fig. U where the left column is three-dimension 
plots and the right column is contour plots. As shown 
in Fig. |l|(a), the DSC order parameter vanishes at the 
core center (dark-blue site) and approaches its zero-field 
value, which is about Ao = 0.08 for the chosen parameter 
values, away from the core center. Notice that the DSC 
order parameter is not uniform beyond the distance £o 
away from the core center. Instead, it is weakly modu- 
lated. The maximum modulation amplitude is less than 
0.05Ao. This modulation is closely related to the ap- 
pearance of the field-induced SDW, as will be discussed 
immediately. Fig. |l|(b) displays the spatial distribution 
of the staggered magnetization of the local SDW order 
defined as M s = (— l) l (n^ — njj_). Clearly, the maxi- 
mum strength of M s is pinned at the vortex core center. 
This AF SDW order exhibits a modulation pattern with 
the satellite peaks (yellow spots) and valleys (dark-blue 
spots) regularly spaced throughout the unit cell along the 
Cu-0 bond directions. The modulation pattern of the 
SDW implies a much longer magnetic correlation length 



as compared to £o, which is consistent with the INS mea- 
surements Q . The appearance of the SDW order around 
the vortex core also strongly affects the electron den- 
sity rii = n i<y- A s shown in Fig. 0(c), at the vortex 
core center, where the SDW amplitude reaches the global 
maximum, the electron density is strongly enhanced. In 
addition, the charge density also exhibits regular mod- 
ulation. By comparing the spatial distribution of the 
DSC order parameter, the field-induced SDW as well as 
the associated CDW, one finds that as long as the abso- 
lute amplitude (\M S \) of SDW reaches a local maximum 
(both yellow and dark-blue spots in Fig. 0(b)), the asso- 
ciated CDW also reaches a local maximum (green spots 
in Fig. 0(c)), that is, a local minimum in the hole den- 
sity, while the DSC order parameter has a local mini- 
mum (shallow red spots in Fig. 0(a)). The one to one 
correspondence between the SDW and associated CDW 
orderings at zero field, was discussed in the context of 
stripes O, |l2[ . To make a closer inspection of 

the periodicity of the SDW and associated CDW modu- 
lation, we perform a Fourier transform of these two quan- 
tities. As shown in Fig. 0(a) , the strongest spectral inten- 
sity of the SDW modulation occurs at the wave vectors 
k = 27r(^;,0) and k = 27r(0, which gives unambigu- 
ously the period of 8ao along the Cu-0 bond directions. 
To verify this point, we have performed a calculation on 
the unit cell of different size N x x N y = 52 x 26, and 
found that the period of the SDW modulation remains 
8ao, indicating convincingly that the periodicity of the 
field-induced SDW modulation is an intrinsic property. 
Therefore, our result explains very well the INS mea- 
sured value of the SDW modulation [Q. We notice that, 
in the optimal doping regime, the circulating supercur- 
rent around the vortex, which favors the isotropy between 
x and y directions, overcomes the one-dimensionality of 
stripes. In the undcrdoping regime, the stripe behavior 
becomes dominant pof . 

Impurity-induced SDW and CDW. The d-wave super- 
conductivity can also be suppressed locally by a strong 
nonmagnetic impurity. The effects of impurities on su- 
perconductors have been of theoretical and experimental 
interest even in its own right for a long time. In view 
of the recent observations of the AF magnetism around 
the vortex core, it is important to see whether a similar 
magnetic structure can also be induced around an impu- 
rity. We perform the numerical calculation on a unit cell 
of N x x N y — 36 x 36 sites. The single-site potential is 
taken to be eo = 100 at the unitary impurity and zero on 
other sites. Since ifiij — in this case, the conventional 
Bloch theorem is used. The results on the spatial distri- 
bution of the three orderings is displayed in Fig. 0. The 
DSC order parameter is depressed dramatically at the 
impurity site and approaches the bulk value at the scale 
£o- Out of this range, almost no modulation of the DSC 
order parameter is seen. The induced staggered moment 
of the SDW is zero at the impurity site and has max- 
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FIG. 1: The three-dimensional (left column) and contour 
(right column) display of the amplitude distribution of the 
c?-wave SC order parameter |A<j| (a), the staggered magneti- 
zation M s (b), and the electron density Sn — ^ CT Ui a — ny (c) 
around one vortex located at the center of an area of 24 x 24 
sites. The size of the whole magnetic cell is 48 x 24. Parameter 
values: t' = -0.2, U = 2.5, V = 1, and ny = 0.84. 
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FIG. 2: The Fourier transform of the spatial modulation of 
the spin density M s (a) and the charge density 8n (b) around 
the vortex. Parameter values are the same as in Fig. El 



ima on the four nearest neighbor Cu sites of the impu- 
rity. Away from the impurity, the induced SDW shows 
a modulation very similar to the vortex case, although 
with a much weaker amplitude. However, the electron 
density, which is zero at the impurity, only exhibits a 
Friedel-like oscillation within a limited range around the 
impurity. The Fourier transform, as shown in Fig. H, 
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FIG. 3: The three-dimensional (left column) and contour 
(right column) display of the amplitude distribution of the 
d-wave SC order parameter |A<j| (a), the staggered magne- 
tization M a (b), and the electron density Sn — m — nj (c) 
around a nonmagnetic unitary impurity (eo = 100) located at 
the center of the unit cell of 36 x 36 sites. Parameter values 
are the same as in Fig. |l| 





BSCCO that show no evidence for CDW modulation 
around the impurities. However, SDW appears to be a 
robust feature induced by the impurity, which should be 
observable in neutron scattering experiments. 

Nuclear magnetic resonance (NMR) measurements 
have shown that when a Cu 2+ in the Cu-0 plane is sub- 
stituted by a strong nonmagnetic impurity, such as Zn 2+ , 
an effective magnetic moment can be induced on the Cu 
sites around the impurity site 22 ^5|. More 

recent NMR measurements |2j| show for the first time 
that near optimal doping, the Kondo screening effect 
observed above the superconducting transition temper- 
ature, is strongly reduced in the superconducting state. 
This indicates the stabilization of the magnetic moments. 
Our result of the strongest staggered magnetic moments 
on the four nearest-neighbor Cu sites of the impurity is 
consistent with the above experiment 25 . This leads 



us to speculate that the magnetism around vortices and 
around impurities may share a common origin. 

In conclusion, we have studied the magnetism around 
vortices and nonmagnetic unitary impurities. In the case 
vortices, the experimentally observed 8ao period of the 
SDW modulation is explained for the first time based 
on a microscopic model. The correspondence between 
the SDW and CDW modulations has also been estab- 
lished. Around the unitary impurity, we have also shown 
the existence of the SDW modulation with a period of 
8ao- The existence of such modulation can be tested by 
neutron scattering experiments. 

Acknowledgments: We wish to thank A.V. Bal- 
atsky, J.C. Davis, B. Khaykovich, G. Ortiz, S. Sachdev, 
OS. Ting for useful discussions. This work was sup- 
ported by the US Department of Energy through the Los 
Alamos National Laboratory. 



FIG. 4: The Fourier transform of the spatial modulation of 



the spin density M a (a) and the charge density Sn (b), around 
a nonmagnetic unitary impurity. Parameter values are the 
same as in Fig. fel 



gives a clear picture: The main modulation period of the 
SDW is 8ao; the CDW is modulated with a period of 2a , 
which is roughly the Fermi wave length. Therefore, the 
correspondence between the SDW and CDW is absent in 
this impurity case. The reason lies in the fact that the 
CDW is very sensitive to the strong impurity scattering. 
This result is consistent with the STM measurements in 
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